A CHAIN RULE FORMULA IN BV AND APPLICATIONS TO 
CONSERVATION LAWS 
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Abstract. In this paper we prove a new chain rule formula for the distributional derivative 
of the composite function v(x) = B(x,u(x)), where u \]a, 6[— ¥ M. d has bounded variation, B(x, ■) is 
continuously differentiable and B(-,u) has bounded variation. We propose an application of this 
formula in order to deal in an intrinsic way with the discontinuous flux appearing in conservation 
laws in one space variable. 
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1. Introduction. In 1967, A.I. Vol'pert in [23] (see also [24]), in view of ap- 
plications in the study of quasilinear hyperbolic equations, established a chain rule 
formula for distributional derivatives of the composite function v{x) = B(u(x)) , where 
u : ft — > K has bounded variation in the open subset fi of M. N and B : R — > R is con- 
tinuously differentiable. He proved that v has bounded variation and its distributional 
derivative Dv (which is a Radon measure on f2) admits an explicit representation in 
terms of the gradient VB and of the distributional derivative Du . More precisely, 
the following identity holds in the sense of measures: 

Dv = VB{u)Vu C N +VB(u)D c u+[B(u + ) - B{u~)\ v u H N ~\j u , (1.1) 

where 

Du = Vu C N + D c u + v u H N ~ l [_J u (1.2) 

is the usual decomposition of Du in its absolutely continuous part Vu with respect 
to the Lebesgue measure C N , its Cantor part D c u and its jumping part, which is 
represented by the restriction of the (N — l)-dimensional Hausdorff measure to the 
jump set J u . Moreover, v u denotes the measure theoretical unit normal to J u , u is 
the approximate limit and u + , u~~ are the approximate limits from both sides of J u . 

The validity of (1.1) is stated also in the vectorial case (see [2] and Theorem 3.96 
in [3]), namely if u : ft — > M. d has bounded variation and B : R d — > R is continuously 
differentiable, then the terms in (1.1) should be interpreted in the following sense: 

Dv = VB{u) -VuC N + VB{u) ■ D c u + [B(u+) - B(u~)] (8 v u [J u . (1.3) 

The situation is significantly more complicated if B is only a Lipschitz continuous 
function. In this case, the general chain rule is false, while a weaker form of the 
formula was proved by Ambrosio and Dal Maso in [2] (see also [21]). 

On the other hand, in some recent papers a remarkable effort is devoted to es- 
tablish chain rule formulas with an explicit dependence on the space variable x . 
This amounts to describe the distributional derivative of the composite function 
v{x) = B(x, u(x)), where B(x, •) is continuously differentiable and, for every s € M. d , 
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B(-,s) and u are functions with low regularity (which will be specified later). These 
formulas have applications, for example, in the study of the L 1 lower semicontinuity 
of approximating linear integrals of convex non-autonomous functional (see [11], [12] 
and [1]). 

The first formula of this type is established in [13] for functions u £ R d ) 
by assuming that, for every s £ R d , B(-,s) is an L 1 function whose distributional 
divergence belongs to L 1 (in particular it holds if B(-,s) £ W 1 ' 1 ^; R d ) ). 

In [12] the formula is proved by assuming that, for every s £ Mr, B(-,s) is an 
L 1 function whose distributional divergence is a Radon measure with bounded total 
variation and u £ W 1 ' 1 ^; R) . 

The case of a function u £ BV(Q) is studied in the papers [11] and [12]. In the 
first paper the authors have established the validity of the chain rule by requiring 
that B(-, s) is differentiable in the weak sense for every set. In the second one it is 
assumed only a BV dependence of B with respect to the variable x . 

The main difficulty of these results consists in giving sense to the different terms 
of the formula. Notice that the new term of derivation with respect to x needs a 
particular attention. For instance in [12] this term is described by a Fubini's type 
inversion of integration order. 

The aim of this paper is to establish a chain rule formula for the distributional 
derivative of the composite function v(x) — B(x, u(x)) , where u :]a, b[— >• K d has 
bounded variation, B(x, •) is continuously differentiable and B(-,s) has bounded vari- 
ation. We assume that there exists a countable set AT C]a, b[ such that the jump set 
Jb(-,s) OI B(-,s) is contained in N for every s e R d . Moreover we require that there 
exists a positive finite Cantor measure A on ]a, b[ such that (D^B)(-, s) <C A for every 
s e K d . For every s £ M. d let ip(-,s) denote the Radon-Nikodym derivative of the 
measure (D^B)(-,s) with respect to A, i.e. 



We show that (see Theorem 4.1 below), under suitable additional assumptions, the 
composite function v(x) := B (x, u(x)) belongs to BV(]a, b[) and for any <j> £ Cq (]a, b[) 
we have 



where u(x + ), u{x-) and B(x + ,s), B(x + ,s) are respectively the right and left limits 
of u and B(-, s) at x . 

The proof is based on a regularization argument via convolutions and on the 
Ambrosio-Dal Maso derivation formula (see [2]). In order to prove the convergence 
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(1.4) 



4>(x)[B(x+,u(x+))- B(x-,u(x-))] , 
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of the regularized terms we follow the arguments as in [12], with the exception of the 
term of derivation with respect to x, which requires a different nontrivial analysis due 
to the possible interaction of the jump points of u and the jump points of B(-, s). 

In order to understand this effect, we consider firstly a piecewice constant function 
u, and we show that, in this case, the contributions of the jump parts can be collected 
as in the summation in (1.4). The general case can be obtained by using a precise 
approximation result, proven in Section 3, of a BV function by piecewise constant 
functions which holds only for functions defined on an interval. By the way, we 
remark that this is one of the technical point where it is crucial the restriction to a 
one dimensional space variable. 

In Section 5, we consider the case d = 1 and we compare our chain rule with 
the formula proven in [12]. We verify the (necessary!) coincidence of the terms of 
derivation with respect to x in the case of piecewise constant functions u. Anyway, 
we remark that the form (1.4) is new also in this one-dimensional case. 

Finally, in Section 6 we discuss the use of our chain rule formula to conservation 
laws with a discontinuous flux. The case of discontinuous fluxes has been intensively 
studied in the last few years (see e.g. [4, 5, 6, 7, 8, 9, 16, 18, 19, 20, 17, 22] and the 
references therein) due to a large class of applications in physical and traffic models. 

We do not address directly the issue of existence or uniqueness of solutions, for 
which wc refer to the references listed above. We remark that the existence results 
are proved only for very special fluxes (tipically, only one jump in the space variable is 
allowed). For what concerns uniqueness, we recall a fairly general result by Audusse 
and Perthame [4], which is based on an extension of the classical Kruzkov method. 

In this framework, using our chain rule formula, we propose a definition of entropic 
solution which is a generalization of the classical one valid for smooth fluxes (see e.g. 
[10]). 

We show that our definition is equivalent, under suitable assumptions, to the 
notion of Kruzkov-type entropic solution obtained using the adapted entropies intro- 
duced by Audusse and Perthame in [4]. Our formula provides a neat environment 
for the treatment of all terms containing a derivative of the composition with a BV 
function which are present in equations of this type. 

We are inclined to believe that the methods here introduced can be useful to treat 
analogous problems in the same context. 

Acknowledgements. The authors would like to thank Gianni Dal Maso and 
Nicola Fusco for stimulating discussions and suggestions during the preparation of the 
manuscript. 

2. BV functions of one variable. In this section we introduce the BV func- 
tions of one variable and we recall the definitions and the basic results (see the book 
[3] for a general survey on this subject). 

We recall that a function u = (u 1 , . . . , u d ) g i 1 (]a, &[; R d ) belongs to the space 
BVQa,b[;R d ) if and only if 

d f b 

TV(u) := sup | / « i D^ i da;:0eC 1 aa,6[;R <, ),||^|| oo <lj<+oo (2.1) 

i=l ^ a 

(if d = 1 the usual notation is BVQa, b[) ). This implies that the distributional deriva- 
tive Du = (Du 1 , . . . ,Du d ) is a bounded Radon measure in }a,b[ and the following 
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integration by parts formula holds: 

r-b r-b 

I u l D<j) l dx = - <pdDu l \/<f>eC 1 {]a,b[;R d ), i = l,...,d. (2.2) 

J a J a 

A measure /x is absolutely continuous with respect to a positive measure A (fj, <C A in 
symbols) if n(B) = for every measurable set B such that X(B) = . We will often 
consider the Lebesgue decomposition 

Du = Vudx + D s u, (2.3) 

where Vu denotes the density of the absolutely continuous part of Du with respect 
to the Lebesgue measure on ]a, b[ , while D s u is its singular part. 

For every function u € BVQa, b[; R d ) the following left and right limits 

1 f x 1 /" x+£ 

u(x-) := lim - / u(y) dy , u(x+) := lim - / u(y) dy (2.4) 

e^0+ e J x _ e e^0+ e J x 

exist at every point x e]a, b[ . In fact, u(x-) is well defined also in x = b, while u(x + ) 
exists also in x — a. The left and right limits just defined coincide a.e. with u and 
are left and right continuous, respectively. 

It is well known that the jump set of u, defined by 

Ju ■= {x e]a, b[: u(x-) ^ u(x + )} 

is at most countable. The singular part D s u of the measure Du can be splittcd into 
the sum of a measure concentrated on J u and a measure D c u, called the Cantor part 
of Du, as in the following formula: 

D s u = D c u + (u(x+) - u(x-))H [J u , (2.5) 

where H° stands for the counting measure. Moreover, we consider the so-called diffuse 
part of the measure Du concentrated on C u -=]a, b[\J u and defined by 

Du := Vudx + D c u, (2.6) 

while 

D*u := (u(.t+) - u(x_))H°lJ u (2.7) 

is called the atomic part of Du. Analogously, we said that a nonnegative Borel 
measure it is a Cantor measure if /j, is a diffuse measure orthogonal to the Lebesgue 
measure. 

If |Du| denotes the total variation measure of Du, we have that |Du|(]a, b[) equals 
the value of the supremum in (2.1); moreover, for every Borel subset B of ]a, b[, 

\Du\(B)= j |Vu|(a;)dx + |D c u|(D)+ ^ |u(x+) - u(z_)| . (2.8) 

Now we recall the classical definition for BV functions of one variable, by means 
of the pointwise variation] for every function u :]a, 6[— >• R d , it is defined by 

n-l 

P V(u) := sup { l u (*i+i) - u ^)l : a < t! < • • • < t„ < (2.9) 
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We remark that every function u having finite pointwise variation belongs to the space 
L°°Qa, b[; since its oscillation is controlled by pV(u). Moreover every bounded 
monotone real valued function has finite pointwise variation and any (real valued) 
function having finite pointwise variation can be splitted into the difference of two 
monotone functions. 

In order to avoid that u changes if it is modified even at a single point, we 
introduced the following definition of essential variation 

eV(u) := sup jpV(v) : v = u a.e. in }a, 6[ j . (2.10) 

Finally, by Theorem 3.27 in [3], the essential variation eV(u) coincides with the 
variation V(u), defined in (2.1). Any function u in the equivalence class of u (that is 
u = u a.e.) such that pV(u) = eV(u) = TV(u) is called a good representative. By 
Theorem 3.28 in [3], we have that u is a good representative if and only if for every 
x e]a, b[ 

u(x) G {(9u(x_) + (1 - 0)u(a+) : 9 e [0, 1]} . (2.11) 

In particular, if (2.11) holds with 9 = (resp. 9 = 1) for every x e]a, b[, we have 
that u = u + (resp. u = u + ), while for 9 = 1/2 u coincides with the so-called precise 
representative 

u'Qr) ;= »(»+) + »(»-) . (2 . 12) 

Any good representative u is continuous in ]a, b[\ J u , and it has a jump discontinuity 
at any point of J u satisfying u(x-) = u(x_), u(x+) = u(x + ). Finally, any good 
representative u is a.e. differentiable in ]a, b[ and its derivative Vu coincides with the 
density of Du with respect to the Lebesgue measure. If not otherwise stated, in this 
paper we always consider good representatives of BV functions. 

For every scalar BV function u the following coarea formula holds (see [14] , The- 
orem 4.5.9): 

rb r+oo r 

/ g(x)d\Du\(x) = / dt g(x)dn°(x) (2.13) 

for every Borcl function g :]a, b[— > [0, +oo[. 

We remark that a Leibnitz rule formula in BV(]a, b[) holds: if v,w £ BVQa, b[), 
then vw € BVQa, b[) and 

D(vw) = v*Dw + w*Dv, (2.14) 

in the sense of measures (see Example 3.97 in [3] and Remark 3.3 in [12]). 

Now we recall the properties of the convolution of a BV function. Let (p be a 
standard convolution kernel and let (tp £ ) E> o be a family of mollifiers, i.e. <p e (x) := 
e~ 1 (p(x/s). For every function u € BVQa, b[; K d ) we define 

u e (x) := (u*ip £ )(x) = / ip £ (x-y)u(y)dy 



for x £]a' ,b'[CC]a,b[ and < e < min(6 — b',a' — a). We have that the mollified 
functions u * tp e converge a.e. to u in ]a, b[ and everywhere in [a, b[ to the precise 
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representative u* (see Proposition 3.64(b) and Corollary 3.80 in [3]). Moreover Vu e = 
V(u * ip £ ) — (Du) * tp £ (see Proposition 3.2 in [3]), where for a Radon measure fi, the 
convolution \i * ip £ is defined as 



{fj,*ip e )(x) := j ip £ (x - y) dfj,(y) . 



Finally, we recall that the measures Vu £ dx locally weakly* converge in ]a,6[ to the 
measure -Du, i.e. for every (j) 6 Cb(]a, b[) we have 



rb pb 

/ (jNu £ dx / <pdDu, as e 

J a J a 



(see Theorem 2.2 in [3]). 

3. An approximation result. In this section we exhibit an explicit piecewise 
constant approximation of a BV function, which is taylored to our needs in the proof 
of Theorem 4.1 . 

Lemma 3.1. Let v <E BV(]a,b[), let J denote its jump set, and let P c]a, b[\J 
be a countable set. Then, for every e > and every finite set P e C P there exists a 
piecewise constant function v e : ]a, b[— > K such that: 

(i) the (finite) jump set J £ of v £ contains all jumps of v of size greater than e / 3 ; 

(ii) TV(v £ ) < TV(v); 

(Hi) J £ n P = and v £ (x) = v(x) for every x G P £ ; 

(iv) v £ (x + ) = v(x + ), v £ (x_) — v(x^), for every x € J f] J £ ; 

(v) \v £ (x) — v(x)\ < e for every x £]a,b[\J (the inequality holds everywhere if v 
is a good representative). 

Proof. Without loss of generality we can assume that v is a good representative. 
Let J — {x^} be the jump set of v. Since v € BV, there exists N e N such that 

5>(4)-«coi<|- (3.i) 

j>N 

Let us define the functions vb,vs'- ]a, b[— > IR by 

( £ [v(xi)-v(x j _)}, iix^{x\...,x N }, 
v B (x) := < xj<x,j<N 

I VsixL) + v(x j ) — v(x J _), if x = x^ for some j < N, 

( H x +) -v(xL)], if x {x l : i > N}, 

v s (x) := < xi<x,j>N 

I vs(x 3 _) + v{x^) — v(x 3 _), if x = x^ for some j > N. 

It is clear from the definition that the functions vb and vs take into account the big 
and the small jumps of v respectively, and that the function vc ■— v — vb — «s is 
continuous in ]a,b[. In addition, vc is uniformly continuous in ]a, b[, since it can be 
continuously extended to [a, b]. Then there exists 5 > such that 

\v c (x) - v c (y)\ < | Vx,y e]o,6[, \x-y\ < S. (3.2) 

Moreover, from (3.1) we have that 

\vs(x)\ < Hx j +)-v(x j _)\ < |, Vxe}a,b{. (3.3) 

j>N 
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Let J £ = {y j }™ =0 , with a = y° < y 1 < ■ ■ ■ < y m = b, be a partition of [a, b] 
satisfying the following properties: 

(a) y l — y l ~ x < S for every i € {1, . . . , to}; 

(b) x j e J e for every jejl,..., iV}; 

(c) for every i G {l,...,m — 1}, if y 4 = for some j G {1, ...,iV}, then 
y l -\y t+1 <£ {x!,...,x N }; 

(d) P n J £ = 0; moreover, each interval ]2/ l_1 ,2/*[ contains at most one point of 
P £ and, in that case, y' ,y % {xi, . . . ,xn}- 

Let i G {1, . . • , to}. For every x, y , j/ 4 [ we have that 

\v(x) - v(y)\ < \v(x) - v c {x) - v B (x)\ + \v(y) - v c (y) - v B {y)\ 
+ \vc(x) - v c {y)\ + \v B (x) - v B (y) \ ■ 

Since ]y l ~ 1 ,y i [ does not contain points of {a; 1 , . . . ,1^} we have that v B (x) = v B (y). 
Moreover, v — vc — v B = vs, hence by (3.2) and (3.3) we obtain 

\v(x) - v(y)\ < \v s {x)\ + \v s {y)\ + \v c (x) - v c (y)\ < e (3.4) 

(x,ye]y i - 1 ,y i [,i = l,...,m). 

Finally, let us define the function v £ : ]a, b[— > R by v e {y % ) = v(y l ) for every i € 
{0, . . . , to}, and, on every interval ]y t ^ 1 , y l [ (i € {1, . . . , to}) by 

(v(x), if0^P £ n]y 4 - 1 ,y 4 hM, 
v £ (x) := I v{y{_), if y l £ {x 1 ,... 7 x N }, 
[^(j/^ 1 ), otherwise. 

It is clear from the construction that (i)-(iv) hold. Moreover, on every interval 
]y 1 ~ 1 ,y l [ (i € {l,...,m}) we have that v e (y l _) = v{y l _) or v e (y\) = v(y\) or 
v e (x) = v(x) for some x &\y t ~ 1 ,y l {, hence from (3.4) we conclude that also (v) holds. 
□ 

Lemma 3.2. Let u G BV(]a, b[; K. d ) be a good representative, let J denote its jump 
set, and let P c]a,b[\J be a countable set. Then there exists a sequence of piecewise 
constant functions u„ e BV(]a, b[, M. d ), neN, satisfying the following properties: 

(i) the (finite) jump set J n of u n does not contain points of P and contains all 
jumps x e J such that |u(a;+) — u(x_)| > 1/n; 

(ii) TV(u n ) < TV(u); 

(in) |u„(x) — u(x)| < C/n for every x €]a, b[ and n G N, where C = 3^/d; 
(iv) for every x £ P there exists n x £ N such that u n (x) = u(x) for every n > n x ; 
(v) for every x G J there exists n x G N such that u n (x + ) — u(x + ), u n (x_) = 

u(x_) for every n>n x . 
Proof. Let P — {zj}j. For every n G N let us apply Lemma 3.1 to each component 
u l , i — l,...,d with e = 3/n and P £ = {z\, . . . , z n }. The conclusion follows from the 
fact that J = \jf =1 J u i and J n = \J d =1 J u i n , n G N. □ 

4. A chain rule formula in BVQa, b{; R d ). Let B :]a, b[xR d ->• R be a function 
such that B (•, w) G BVQa, b[) for all w G R d . We recall that for every w G R d 

(D x B)(-,w) = (V x B)(.,w)dx + (D c x B)(.,w)+ ^ [B(x+,w) - B(x-,w)]5 x (4.1) 

is the usual decomposition of the measure (D x B)(-,w) with respect to the Lebesgue 
measure, where A/" w := J_b(-.w) lS the jump set of £>(•, w) . 

Theorem 4.1. Let B :]a, b[xR d — >• R &e a locally bounded function such that 
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(Al) for all w e R d the function B(-,w) belongs to BV(]a,b[) and there exists a 
countable set N c]a, b[ such that for every w e M. d we have 

A/w C TV; 

f-Ai?,) for every compact set M CM. d there exists a finite positive Borel measure fiM 
in ]a, b[ such that for every w, w' e M and every Borel set A C]a, b[ 

\{D x B){-,w) - (D X B)(-,W)\(A) < \w-w'\ m (A); 

(A3) for all x £]a,b[\J\f the function B(x,-) belongs to C 1 (R d ) and, for every 
compact set M C R d , there exists a constant D M > such that 

\(D w B)(x,w)\<D M , Vxe]a,b[\N,w eM; 

(A4) the function (D^B) (•, w) belongs to BV(]a,b[;R d ) for every w e M. d ; 
(A5) there exists a positive finite Cantor measure X on]a, b[ such that (D X B)(-, w) <C 
A for every w e M. d . 

Then for every u e BV(]a, b[; R d ) the composite function v(x) := B(x,u(x)), x e 
]a,b[, belongs to BV(]a,b[) and for any </>e C^Qa, b[) we have 



/ 4>'{x)v{x) dx = — / (j){x)(y x B){x, u(x)) dx 

J]a,b[ J]a,b[ 

(j)(x)ip(x, u(x)) dX 



a,b[ 



/ 4>{x)(D vi B)(x,u(x)) ■ \7u(x)dx (4.2) 



<f>(x)(D w B)(x,u(x)) -dD c u(x) 

a,b[ 

- <Kx)[B(x + ,u(*+))-B(x-M*-))], 

xeA/UJ u 

where for every w e R d the function ip(-,w) is the Radon-Nikodym derivative of the 
measure (D x B)(-,w) with respect to X, i.e. 

Remark 4.2. By (A2) we obtain that for every compact set M C R d there exists 
a constant Cm such that 

\(D x B)(-,w)\(]a,b[)<C M , VweM. (4.3) 

Moreover, for a.e. x G]a,b[ we have that 

\(V x B)(x,w)\<C M , VweM, (4.4) 

and, for every x e]a, b[, 

\ip{x,w)\ < C M , VweM. (4.5) 
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In addition, for a.e. x E]a,b[ and for every Borel set A C]a, b[ the functions w i->- 
(V x B)(x, w) and w h- > (D x B)(-,w)(A) are Lipschitz continuous in K. d . Finally, for 
X-a.e. x e]a, b[ the function w i->- ij)(x,Vf) is Lipschitz continuous too. 

Remark 4.3. Formula (4.2) can be rewritten in a more explicit way as 



/ (f>'(x)v(x) dx = 

J]a,b[ 

- (f)(x)(V x B)(x,u(x))dx - / <p(x)(D vr B)(x,u(x)) ■ Wu(x)dx 

J]a,b{ J]a,b{ 



- (f>(x)i{>(x,u(x))d\ - / (f)(x)(D vr B)(x,u(x)) ■ dD c u(x) 

J]a,b[ J]a,b[ 

_ Sr- Mx) \ B(x + ,u(x+)) + Bjx+Mx-)) _ B(x-,u(x + ))+B(x-,u{x-)) 

- £0(a)[fl>,u(a^)) u(a;_))] , 



where for every x G]a, b[ and w &M. d 

B * [x , w) ;= g(» + ,w)+g(»-,w) 

is the precise representative of the BV function x h-» B(x,w) . 
In fact, it is easy to check that for every x e J U C\JV we have 

B{x+, uO+)) - B(x-,u(x-)) = [B*(x, u(x+)) - B*{x, u(z_))] 
B(x + ,u(x + )) + B(x + ,u(x-)) B(x-,u(x + )) + B(x-,u(x-)) 
2 2 

in particular, for every x e J u \Af we have 

B(x+,u(x+)) - B(x-,u(x-)) = B(x, u(x+)) - B(x, u(z_)) 

= B*(x,u(x+))-B*(x,u(x-)), 

and for every x € N \ J u we have 

B(x + ,u(x + j) — B(x-, u(x_)) = B(x +1 u{x)) — B(x-,u(x)) 
_ B(x + ,u(x + )) + B(x + ,u(x-)) B(x-,u(x + )) + B(x_, u(x_)) 



(4.6) 



Proof of Theorem 4.1 

Since the proof of Theorem 4.1 is rather long, it will be convenient to divide it 
into several steps. 

In Step 1, following the regularization argument of Ambrosio-Dal Maso (see [2]), 
we consider the mollification B E (x,w) of B(x,w) with respect to the first variable. 
We observe that, for every test function <f> e C^Qa, b[), the integral 

/ 4>'(x)B e (x,u(x))dx 

J]a,b[ 
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converges to the left-hand side of (1.4) as e — > + (see (4.8)). Then, for e small 
enough, we decompose this integral (using the chain rule formula for C 1 functions) as 



/ (j)'(x)B s {x, u(x)) dx = / <f>(x) D w B e (x, u(x)) dDu(x) 

J]a,b[ J]a,b[ 

+ <f>( x ) [B s {x,u(x+)) - B e (x,u{x-))} 



f 4>(x) (D x B e )(x, u(x)) dx =: D s + J s + I e 

J]a,b[ 



and we study the convergence of each one of the three terms D s , J s , I e appearing at 
the right-hand side as e — > + . 

The limits of D e and J e are computed respectively in Steps 2 and 3 following the 
lines of [12]. 

The limit of I e is far more difficult to analyze, because of the possible interaction 
between the jump set of u and the jump set of B(-,u). In Step 4 we compute this 
limit in the special case of u piecewise constant. Finally, the general case is proved 
in Step 5 relying on a carefully chosen approximation of a BV function by means of 
piecewise constant functions, whose construction has been shown in Lemma 3.2. 

Step 1. Fix 4> e C*(]a, b[) and let ip e — ip e (x) be a standard family of mollffiers. 
Let us define 



B £ (x, w) := / ip £ (x - y) B (y, w) dy 

J]a,b[ 



for x e]a', b'[ and w e R d , where supp <fi C [a', b'] c]a, b[ , and < e < min{b — b',a — 

<<'!• 

We claim that B e e C' 1 (]a' , b' [xR d ) . Firstly we prove that D X B £ is locally Lip- 
schitz continuous in }a',b'[xR d . In fact, by hypothesis (A3) for every compact set 
D C]a' ,b'[xR d and for every (xi,wi), (x2,w 2 ) G D there exists a constant Cd such 
that 



(AAXa^Wi) - {D x B £ )(x 2 ,w 2 )\ 

= / [<P' E { X 1 ~y)B{v^i) -V's{x2 -y)B(y,w 2 )] dy 

J]a,b[ 

/ [<P , E { x i-y)-<Pe( x 2-y)]B{y,w 1 )dy 

J]a,b[ 

+ / ip' £ {x 2 -y)[B(y 7 w 1 )- B(y 7 w 2 )]dy 

J]a,b[ 



< 



< —(\xi - x 2 \ + |wi - w 2 |) . 



Moreover, we prove that D^B,, is continuous in ]a', b'[xR d . In fact, for every sequence 
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(x n ,w n ) converging to (x, w) in }a',b'[xW i we have 
|(D W -B e )(a: n ,w n ) - (D w B £ )(x, w)| 

= / [<Ps(x n - y)(D w B)(y,Mv n ) - tp s (x - y)(D Vf B)(y,w)] dy 

J]a,b[ 

< / \<p e (x n - y) - (fe(x - y)\ \(D vf B)(y, w„)| dy 

J]a,b[ 

+ / \<Ps{x - y)\ \(D w B)(y,w n ) - (D Vf B)(y,w)\ dy . 

J]a,b[ 

The first integral tends to 0, as n — > oo, since by (A3) 

f C 

/ \<Pe{x n - y) - ip e (x - y) \ \(D w B)(y,w n )\ dy < —D M \x n - x\ , 

J]a,b[ e 

and the second one tends to 0, as n — > oo, by the continuity of the function {D vl B){y 1 ■) 
for a.e. y £]a, b[, the boundedness of D W B and by the Lebesgue dominated conver- 
gence theorem. 

Let u e BV(}a 7 b[;M. d ) and define 

v £ (x) := B £ (x,u(x)), xe]a',b'[. 

Since B £ e C 1 (]a', b'[xR d ) we can apply the chain rule formula (see Theorem 3.96 in 
[3]) to the composition of the function B £ with the BV map x (x, u(x)), concluding 
that v £ e BV(]a',b'[) and 

/ (p'(x)v £ (x) dx = - 4>{x) (D x B £ )(x,u(x)) dx 

J]a',b'[ J\a',b'\ 

- / (j){x)(D u B £ )(x,u{x)) -dDu(x) 

J]a',b'[ 

- J2 0(x)(B £ (x,u(x + )) - B £ (x,u(x-))) 
xe j u n]o',6'[ 

= -/ ((>(x)(D x B £ )(x,u(x)) dx (4 ' 7) 

J]a',b'{ 
d 

-J2 4>(x)(D wt B £ )(x,u(x))dDu*(x) 
i=iJw,V[ 

-E E Hx)(u l (x + )-u\x-)) f (D wi B £ )(x,w s (x))) da, 
*=i xeJ,An]a' b'l 



»n]c 

where Du and flu' denote the diffuse parts of the measures Du and Du l respectively, 
and w s (x) := u(x_) + s(u(a; + ) — u(a;_)). 

Since B is locally bounded and the functions B £ (-,w) converge a.e. in ]a',b'[ to 
B(-,w), by Lebesgue dominated convergence theorem we get 

lim / (j>'(x)B £ (x,u(x))dx= <j>'(x)B (x, u(x)) dx . (4.8) 

e ^ 0+ J]a,6[ J]a,b[ 
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Step 2. We shall prove the convergence of the diffuse part, i.e. for every i = l,...,d 
we prove that 

lim / ^{D^B^x^x^dDv? = [ cf)(x)(D w iB)(x,u(x))dDu\ (4.9) 

£ ^ 0+ J]a'.b'[ J]a,b[ 

Using the coarea formula (2.13), we get 

/ <t>{x)(D w iB e )(x,u 1 {x),...,u i {x),...,u d (x))dDu i (4.10) 

J]a',b'[ 

= / c/ ) (x)(D w ,B £ )(x,u 1 (x),...,u l (x),...,u d (x))—-(x)d\Du l \ 
J]a',b'[nc ui \Du z \ 

r+oo r f) A 

= dt <f>(x) (D w i B e ){x, u(x))——r- (x) dH° 

J -co JR<t<«;}nC B , \ Uu I 

dt / <j>{x){D wi B e ){x,u\x), u d {x))——{x) dH° . 

-oo J{u'=t}nc ui \ Du \ 

Now, by (A4) we have that for every i = 1, . . . , d and for every w e M d 

(D Wt B e )(x,w) -> (D Wi B)*(x,w) Vx e]o,6[ (4.11) 

as £ — >• 0. Therefore, for a.e. i S R, we have 

lim / 0(a:)(D 1u «B e )(a:,« 1 ,...,t,...,u <, ) T77ir dW 

= / <l>{x){D w iB)*(x, u 1 , . . . ,t, . . . , u )— -r- cfH° . 

From this equation, using the local boundedness of (D Wi B)* and the fact that, by the 
coarea formula (2.13), 

/+oo 
H° (K = t} n C ui ) dt = \Du l \(C u >) < oo , 
-oo 

we can pass to the limit in (4.10) and by Lebesgue dominated convergence theorem 
we get 

lim / (j)(x)(D wi B £ )(x,u 1 {x),...,u i (x),...,u d {x))dDu i 

J]a'.b'[ 

/+oo ,■ 
dt / 4>{x)(D u iB)*(x, u 1 , . . . ,t, . . . , u d ) dDu 1 . 

-oo J{ M i=t}nc ui 

From this equation, using the coarea formula (2.13) again, we immediately get (4.9). 

Step 3. We shall prove the convergence of the jump part, i.e. for every i = l,...,d 
we prove that 



lim + ( t>(x)(u l (x + ) - u*(x_)) [ (D w iB E )(x,w s (x)))ds 

= Yl Hx)(u\x + )-u\x-)) j (D wi B)*(x,w s (x)))ds, 
xeJ ui r\]a',b'{ 



(4.12) 
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where w s (x) := u(x_) + s(u(x + ) — u(x_)). Let us fix i £ {1, . . . , d}, and let J' ui := 
J u iC\}a', b'[— {j/ J '}jeN- F° r every h £N there exists fc(/i) <E N such that 



j=k(h)+i 

Then the following estimate holds: 



oo 

£ | u v + )-«Vi)l<^ 



/.l 

GJ', ^° 

< IMUE !"'(»+) " u Vj| / (^Be)(2/ j ,w s (^))-(£> wiJ B)*(^,w s (^)) 

+ W=o E |^V+)-« < (yi)| / |(^Be)(^,W s (^))-(^B)*(^,W s (^)) 

< Halloo / y"|« i (l4-)-«Vj|(^5 e )(^\w a (^'))-(£) TO iB)*(^,w a (^)) 

Jo - n 



+ 2C||^|| 00 2 lu'^J-uVJl, 
i>fe(h) 

where C := \\D w iB\\ Loo Q al tb ,^_ M ^ and ||u 4 ||oo < M. By Lebesgue dominated 
convergence theorem, the first integral is infinitesimal as e — > 0, since (D w iB e )(x,w) 
and (D w iB)*(x, w) are locally bounded functions and for every x £]a' , b'[ and w G R d 
we have that (D w iB e )(x,w) — > (D w i £?) * (x, w), as e — > 0. Therefore, letting first e 
tend to zero and then h tend to oo, we immediately obtain (4.12). 

Step 4. In this step, we consider a piecewise constant function u: ]a, &[— > R d of the 
form 

N 



"(x)=J2v i X * [aim+l] (x), 



i=0 



where v°, . . . , v N £ R d , a = a < a\ < . . . < ajy < a^+i = b and we prove that 
lim / <j)(x) (D x B e )(x, u(x)) dx 

e ^°J]a'M[ 

= [ <j>(x)(V x B)(x,u(x))dx + [ ^) d(D ^ )( '' u) dA 

J]a',b'[ J]a',b'[ « A 

B(x + ,u(x + )) + B(x + ,u(x-)) B(x-,u(x + ) + B(x-, u(x_)) 



+ E 

xeJVuj u 



(4.13) 



In order to simplify the notation, let us denote by \% the characteristic function 
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X[o»,o i+ i]- From the very definition of u and (2.14), we have that 
L 



:= / (j)(x) (D x B E )(x,u(x))dx 

J]a',b'[ 
N 

= E/ <j>(x)x;(x)(D x B £ )( X y)dx 

N 

= -V/ B e (x,v l )dD(^) 

N 

= -V/ ( f> , (x) X i(x)B e (x,v i )dx 

N 

+ Y1 [<t>{ai+i) B £ {a l+1 y) - c/)(a l )B e (a i ,v i )] 



i=0 



Passing to the limit as e —¥ we obtain 



N 

lim/ £ =-V/ (j)'(x)xi{x) B(x,v l )dx 



i=0 
N 



(4.14) 



+ Wi+i) B*(a l+l y) - 4>{ ai ) B*(a t ,v*)] 



i=0 



Let us consider the integrals at the right-hand side of (4.14). Using again (2.14) we 
have that 



/ ( f>'(x)xi(x)B(x,v i )dx= [ ^dDixiBi-,^)) 

J]a,b[ J]a,b[ 



a,b[ 



( />(x)x*(x)d(D x B)(;v*) + 



a,b[ 



(x) B*(x,V)dD Xl 



= f c/ ) (x)xHx)d(D x B)(;v l ) + cl ) (a l )B*(a l ,v l )-4>(a l+1 )B*(a l+1 ,v l ). 

J]a.b[ 

Substituting this expression into (4.14) we thus obtain 

N 

I:=limI £ = J2 / Hx)xUx)d(D x B)(-,v l ). 

Finally, let us decompose each measure (D x B)(-,v l ) in the canonical way (4.1). It is 
not difficult to check that 



]T / ftx) X* (x) V x B{x, v') dx + £ f cj>{x) X * (x) dD * B ^ d A 



dX 



N 



+ ]T <Kx)[B{x+,u(x))-B{x-Mz))] 



B(a t +, v J ) + B( ai +, v 4 - 1 ) B(ai-, v l ) + B(oj-, v^ 1 ) 1 



(4.15) 
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The first two terms coincide respectively with 

4>(x) S7 x B(x, u(x)) dx , / (j>(x)tp(x,u(x))d\. 

a,b[ J]a,b[ 

The last two summations take into account the jump points x € J u and x € Af \ J u , 
respectively. Again, it is not difficult to check that, in both cases, the corresponding 
term can always be written as 

B(x + ,u(x + )) + £?(x+,u(x_)) B(x-, u(x+) + B(x-, u(x_)) 



2 2 
so that (4.13) follows. 

Step 5. In this step, we shall prove that formula (4.13) holds for every function 
u e BVQa, b[; R d ), i.e. we prove that 

(4.16) 
(4.17) 



I = lim I s 

6^0 + 



where 



:= / (f>(x) {D x B s )(x, u(x)) dx 

J\a',b'[ 



and 



/:= / (j>(x)(V x B)(x,u(x))dx + / </>(x)ip(x, u(x)) dX 

J]a,b[ J]a,b[ 

+ ~ B ( X+,U ( X+ ^ + B ( X+,U ( X ~^ - B ( x -> u ( x +) + B(x-,u(x-)) 

Let u G BV(]a,b[;R d ) and let (u„)„ be the sequence of approximating piecewise 
constant functions given by Lemma 3.2 with P = Af \ J u . 
Fixed e > 0, we set 



Ie ■= I 4>{x){D x B e )(x,u n {x)) dx. 

J]a'b'\ 



(4.18) 



By Lebesgue dominated convergence theorem and the continuity of (D x B s ){x, ■ 
(which follows by B £ € C^Qa', b'[xR d )), for every e > we have that 



I £ = lim I? . (4.19) 

n— >oo 

More precisely, we claim that 

\Ie - 41 < IK - u|| 00 /i M (]a,6[)||^|| 00 Ve > small enough. (4.20) 
Namely, by hypothesis (A2) we have 



\I?-I e \ < 



[ \(/>(x)\\ f ip e (x-y)d\(D x B)(;u n (x))-(D x B)(;u(x))\(y) 

J]a,b[ l J]a,b{ 



dx 



< \<j>(x)\\u n (x) - u(x)| / ip £ {x - y)djj, M 

J]a,b[ l J]a,b[ 



(v) 



dx 



<IK ~ u lloo||0|| c 



a,b[ 



a,b[ 



ip £ (x-y)dx djjL M {y) 



=|K - ullooll^Hoo/XMQa, b[) . 



(4.21) 
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On the other hand, by Step 4 we have for every n E 



lim r = r 

e->0+ 



(4.22) 



where 



and 



I 71 := <j>{x){V x B){x,u n {x))dx + / <j)(x)tp(x, u n (x)) dX + S n , (4.23) 

J]a,b{ J]a,b[ 

■B(x+,u n (x+)) + B(x+,u n (x-)) 



S n := J2 ^ 



-+ 



B(x-,u n (x + ) + B{x-,u n {x-)) 



We claim that 



By Remark 4.2 we have that 



/ = lim I n . 

n— >oo 



(4.24) 



lim / 4>{x)(S7 x B){x, u n (x)) dx = / (f>(x){y x B)(x, u(x)) dx 

n ^°°J]a,b[ J]a,b[ 



and 



lim / <fi(x)il)(x, u n (x)) dX = / (p(x)tp(x,u(x))dX. 

n ^°°J]a,b[ J]a,b[ 

It remains to show that S n converges, asm +oo, to 

S '= ^ ' B(yX+,U ( X+ ' > ' > + B ( x +' u ( x -^ _ B(x-,u(x+) + B(x-,u(x-)) 

Lct h = U„ J u„ , h = N\h and I = h Ul 2 ■ We recall that P = Af\J u and J u C h ■ 
Since, by construction, J Un n P = for every n £ N, we have that 7i fl P = and 
N U J u C I. Hence both summations in S n and S can be extended to the bigger set 
I = {xi}, since it is easy to check that the added terms are all zero. Thus we can 
write 



where 

cii := <f>{x l ) 



B(4, u n (4)) + u n (z!_)) UnK) + M*-)) 



2 2 
S^^u^V)) +P«,u(a;!_)) B(x!_, 11(24) + B(xt, u(a;L)) 



Let i? > maxjUunlloo, IjuHoo} and let M = Br{Q). From assumption (A2) we have 
that 

|P(z + ,w)-P(a;_,w)| < |P(a; + ,0)-P(a;_,0)|+P/XM({a;}), x €}a,b[, weM. 
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Since 

Kl < h := H^Hoo [\B{x\,Q) - B(xt,0)\ + Rm({xi})] , Vn e N , 

and 

J2 b * < H\\oc{TV(B(;0)) + RuMQ^bl)} , Vn e N , 

in order to prove that lim„ 5„ = 5, by Dominated convergence theorem, it is enough 
to prove that lim„ a™ = dj for every i e N. 

We have three cases. If Xj £ P = Af\ J u , then u and every u„ are continuous at 
Xi. Moreover, for every n large enough, u ra (xj) — u(xj), hence a™ = a^. If Xj G J u , 
then for every n large enough we have that u„(xi+) = u(x{+), u„(xj— ) = u(xj— ), 
hence again a" = a*. Finally, let us consider the case Xj e 7i \ J u . Since Xj A/", the 
function B(xi, •) is continuous in R d . Moreover, since Xj ^ J u , also u is continuous at 
Xi and u„(xj+), u„(xj— ) — > u(x,), so that lim n a™ = aj. Therefore (4.24) is proved. 

In order to prove (4.16), let us fix r\ > 0. By (4.20) and by (4.24) there exists 
n € N such that 

1 7^° - I £ \ < | Ve > small enough 

and 

|7™° -7| < \. 
Moreover by (4.22) there exists £o > such that 

|/no _/no| < | V0<e<e o . 

Then 

\h ~I\< \h - Ie°\ + l 7 e"° ~ + I 7 ™" ~ 7 I < »? V0 < £ < £ . 

Therefore (4.16) is proved and this concludes Step 5. 

Finally, the thesis of the theorem is obtained by collecting all the Steps. □ 

In view to the applications to conservation laws (see Proposition 6.1) we need 
to generalize formula (4.2) in order to integrate a BV function with respect to the 
measure (B(x,u(x))) x . 

Corollary 4.4. Let B :]a, fo[xR — >• R be a function satisfying the same assump- 
tions of Theorem 4-1- Let g :]a, 6[— > R be a BV -function such that J g C J\f . 

Then for every u e BVQa, b[; R d ) and (f) e CgQa, b[) we have 

[ 4>{x)g*(x)d{B(x,u(x))) x = [ 4>(x)g(x)(V x B)(x,u(x))dx 

J]a,b[ J]a,b[ 

+ / cj)(x)g(x)tp(x, u(x)) dX 

J]a,b[ 

+ <p(x)g(x)(D w B)(x,u{x)) ■ Vu(x)dx (4.25) 

J]a,b[ 

+ / <p(x)g(x)(D v ,B)(x,u{x)) -dL» c u(x) 

J]a,6[ 

+ Yl 0(x).9*(x)[i?(x + ,u(x + ))-B(x_,u(x_))] . 
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Proof. Let g £ — g * ip £ be the standard mollified functions of the BF-function g. 
We recall that g e pointwise converges (everywhere) in }a, b[ to the precise representa- 
tive g* , as e — > . We apply Theorem 4.1 by using 4>{x)g e (x) as test function. The 
conclusion follows by Lebesgue dominated convergence Theorem. □ 

In the next corollaries we consider B(x,w) with a particular structure. 
Corollary 4.5. Let K g BVQa,b[) and f g C 1 ^) . T/ien /or every u g 
BVQa, &[; R rf ) tfte function v :]a, &[— > R ; defined by 

v(x) := K{x)f (u(x)) , xg]a,i)[, 

belongs to BV(]a,b[), and for any <f> g CgQa, 6[) we /iat>e 

/ 4>\x)v(x) dx = - [ <f>(x)(f(u))*(x)dDK(x) 

J]a.b[ J]a,b[ 

-[ <f>(x)K(x)(Vf)(u(x)) -Vu(x)dx 



(4.26) 



a,6[ 

0(x)if(x)(V/)(u(x)) -dD c u(a;) 

a,6[ 

- ^0( a ;)A'*( a: )[/(u(x + ))-/(u(a : _))] ) 

w/iere (/(u))* and if* are f/ie precise representatives of the BV functions /(u) and 
if respectively. 

Proof. It is sufficient to observe that the function B(x,w) := K(x)f(w) satis- 
fies all the assumptions of Theorem 4.1 . For instance, hypothesis (A2) is satisfied 
since for every compact set M C M. d we can choose [im '■= Cm\DK\, with Cm = 
maxtugM |V/(w)|, and hypothesis (A5) is satisfied since we can choose A := \DK\ . □ 

Corollary 4.6. Let f : MxR d -> R and K :]a,6[-» R fee two functions satisfying 

(i) K e BV(]a,b[); 

(ii) the function f = f(y,w) belongs to C^R x R d ) ; 

(mj £/ie function x ^ f w (K(x),w) belongs to BVQa, b[) for every w g R d ; 
(iv) for every compact set Dclx R d there exists a constant Ld such that 

\f y (y,w) - f y (y,w')\ <L D \w-w'\ V(y, w), (y, w') eD 

and 

l/w(y,w)-/ w (y',w)| <L D \y-y'\ V(y, w), (t/, w) g £) . 
T/ien /or every u g BV(]a, 6[; R d ) i/ie function v :]a, 6[— > R, defined by 

v(x) := f (K(x),u(x)) , xg]a,6[, 
belongs to BVQa, b[), and for any <fi g Cg(]a,6[) we /lave 

/ 0'(2>(x)a!x = - / 0(a;)/ w (if(a;),u(x)) • dDK(x) 

J]a,b[ J]a,b[ 

-[ Hx)U(K(x),u(x)) -dDu(x) (4.27) 

J]o,6[ 

- ^ 0(x)[/(if(x + ),u(x + ))-/(if(x_),u(.x_))] . 
ie JuUJk 
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Proof. We observe that the function B(x,w) = f(K(x),w) satisfies all the as- 
sumptions of Theorem 4.1 . In particular, AT = Jk, 

D c x B(.,w) = f y (K(x),w)D c K<^ \D C K\=\ 

and V(*,w) = f y (K(x),vr)$%(x). □ 

5. Comparison with other chain rule formulas. In [12] it was proved a 
chain rule formula for function u : M. N — > M. In Theorem 5.1 we recall this formula 
which coincides to formula (4.2) in the case of d = 1 = N. Although the two formulas 
look like very different, we explicitely show that they concide for piecewise constant 
functions. 

Theorem 5.1. Let B :]a, b[xR — > R be a locally bounded B or el function. Assume 
that B(x, 0) = for all x e]a, b[ and 

(i) for allteR the function B (•, t) <= BVQa, b[) ; 
(ii) for all x G]a, b[ the function B (x, •) belongs to C 1 (R) ; 
(Hi) the function DtB is locally bounded, for all t G K the function (DtB) 
belongs to BV(]a,b[) and for every compact set M C M 



/ 



\D x (D t B)(;t)\(]a,b[) dt<+^. 

Then, for every u G BVQa, b[) the composite function v(x) := B(x,u(x)), x £}a,b[, 
belongs to BVi oc (]a,b[) and for any (f) G Cg(]a, b[) we have 

p r+OO r 

/ <j)'(x)v{x)dx=4 dt sgn(t)xn ut (x)4>(x)dD x {D t B)(-,t) (5.1) 

J]a,b[ J-oo J]a,b[ 

f (f>(x)(D t B)(x,u(x))Wu(x)dx 

J]a,b[ 

4 (t>(x)(D t B)(x,u(x))dD c u(x) 

J]a,b[ 

$2<f>(x) [B*(x,u(x+)) — B*(x,u(x-))] , 

xEJ u 

where Uj t = {x G]a, b[: t belongs to the segment of endpoints and u(x)} and Xq u 
and B*(-,t) are, respectively, the precise representatives of the BV functions xo„ t 
and B(-, t). 

Proof. It is a consequence of Theorem 1.1 in [12], with N = 1 and 



B(x, t) = b(x, s) ds . 
Jo 



We recall that in our case for the approximate limits u + (x) and u (x) we have u + (x) = 
v u {x)u(x + ) and u~ (x) = v u (x)u(x-), where v u = ±1 is the normal at a jump point . 
□ 
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Remark 5.2. When d = 1, if we assume also all the hypotheses of Theorem 4-1 
and we compare formulas (5.1) and (4-2), we conclude that 



/ + OC P 
dt se^(t) X ^ t (x)4>(x)dD x {D t B)(;t) = 
-co J]a,b[ 

(j){x)(y x B){x, u(x)) dx + / 4>(x)ij){x, u{x)) dX 

a,b[ J]a,b[ 

B(x+,u(x+)) + B(x + , u(x-)) B(x-,u{x+)) + B(x-,u(x-)) 



(5.2) 



For piecewise constant functions this formula can be proved by using formula (4-15) 
and the following proposition. 

Proposition 5.3. For every piecewise constant function u : }a, £>[—>• M of the form 



N 



U ( X ) = E W< X*a,,a, + i]( a 0. 



i=0 



where v°, . . . , v N G M, a = a < a\ < . . . < ajy < ajy + i = b, we have that 



f °°dtf S gn(t)xh u , t (x)<f ) (x)dD x (D t B)(;t) = f2 [ ^)x\ auai+l] {x) d(D x B)(-, „*) . 

J — oo J]a,b[ j_q J]a,b[ 

Proof. It is not restrictive to assume that u > 0. By the Leibnitz formula (2.14) 
we have that 

/+oo P 
dt / sgn(i) X ^_ t (x)0(x) dD x {D t B){; t) 
-oo J]a,b[ 

/■+OO p p + OQ p 

= dt ct>(x)dD x ( X n u , t ( D tB))(;t)- dt ct>(x)(D t B)*(x,t)dD X n u , t ( x ) 

Jo J\a,b{ JO J]a,b{ 

=: h+h- 

Since B(-,0) — 0, we have that 



h 



r + oo r 

= dt <t>'{x) X n u , t ( x )( D tB)(x,t)dx 

Jo J] a M 

= - [ <f/{x) { [ ( ] (D t B)(x,t) dt] dx 

J]a,b[ \J0 J 

= — 4>'{x) B(x, u(x)) dx 

J]a,b[ 

= ~J2 <P'(x)x] ai , ai+1 [(x)B(x,v l )dx 
l= o ■'KM 



N 

E 

1=0 



/ <l>(x)x* [ai , ai+l] (x) d(D x B)(x, v*) - cj>(a l+1 )B*{a l+1 y) + ^B*^) 

J]a,b[ 

(5.3)' 
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For what concerns the second term I 2 , let us observe that D\Q u t is an atomic measure 
with support contained in {a\, . . . , ajy}. Moreover 



Therefore, by Fubini's theorem we obtain 



1, if Vi-i < Vi and t e [vi-i,Vi[, 
-1, if > Uj and t e 
0, otherwise . 



AT 



12 = - X] 



i=l 



{D t B)*{a u t)dt. 



We claim that 



so that 



/ * (D t B)*(ai,t)dt = B*(ai,Wi)-B*(ai,Wi_i; 

Vi-l 



N 



h =-20(a i )[B*(a < ,i; i )-B*(a < ,i; < _i)]. 



(5.4) 



Namely, since G W(]a,6[), (D t B)(-,t) e W(]a,6[) and |(D(B)(i,t)| < £> M 

for every x e]a, 6[, we have that 



lim (D t B)(x,t)+ lim (D t B)(x,t) 

X— >CLi+ x—tai — 



dt 



/ (D t B)*( ai ,t)dt=- 

J Vi-l ^ JVi-1 

1 r* 1 r* 

= lim - / (D t B)(x,t)dt + lim - / {D t B)(x,t)dt 

= lim lim i[5(x, v l ) - 5(x, v^ 1 )] 

2;— >ai+ 2 >ai— 2 

= B*(a,i,Vi) - B*(a,i,Vi-i) . 

Finally, the conclusion follows from (5.3) and (5.4). □ 

6. An application to conservation laws. In this section we shall apply the 
chain rule formula in order to study a scalar conservation law where the flux depends 
discontinuously on the space variable: 



i t {x, t) + B(x, u(x, t)) x = 0, (x, i)elx [0, +00) 



(6.1) 



where B:lxR^Risa function satisfying the assumptions of Theorem 4.1 (with 
]o,6[=R). 

For every x £ R we define the set of pairs (ii-,u + ) satisfying the Rankine- 
Hugoniot condition 

A x = {(u-,u + ) eMxI: B{x-,u_) = B(x + ,u + )} . 

We define an entropy-flux pair (r],q) associated to (6.1), as a pair of functions 
77, q : M x M -> R such that: 
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(El) for every x G R the function fy(x, •) is convex and r) u is locally bounded in 
R x R; moreover, for every u G R the functions r/(-,w), r) u (-,u) belong to 
BV(M) and their jump set is contained in Jr-; 

(E2) q(x,-) e Lip loc (M) for every x G R, and q(-,u) G W(R), J g (.,„) C J K for 
every it € R; 

(E3) rj u (x, u)B u (x, u) — q u (x, u) for every x € R \ Jr- and u G R; 

(E4) q(x + ,u + ) — q(x-,u-) < for every x G R and every (u-,u+) G A x . 
Proposition 6.1. Let u be a bounded piecewise C 1 solution of (6.1) with J u (-,t) ^= 
Jk for every t G]0,T[ ; and let (r],q) be an entropy- entropy flux pair associated to (6.1). 
Then u satisfies the following inequality 

{V(x,u)) t + (q(x,u)) x <0 (6.2) 

in the sense of measures, i. e. 

[ [ <t>(x, t) d [(rj(x, u(x, t))) t + (q(x, u(x, t))) x ] < (6.3) 

JjRx]0,T[ 

for every function cj> : Rx]0, T[— >• [0, oo[ continuous with compact support. 

Proof. We remark that, by the chain rule formula and since rj does not depend 
on f, we have 

(r)(x,u(x,t))) t = r)Z(x,u)u t {x,t)dt 
in the sense of measure. Then 

// <p(x,t)d(n(x,u)) t = <p(x,t)f]l(x,u)u t (x,t)dt, 

J JRx]0,T[ J JRx]0,T[ 

so that 

// <j)(x,t)d(T)(x,u)) t + (q(x,u)) x 

J Jlx]0,T[ 

= // <t>(x,t)(ri*(x,u))ut(x,t)dt+ (p(x,t)d(q(x,u)) x , 

jiix]o,r[ J Jmx]o,t[ 

where (t] u {x, u))* is the precise representative of the composition of f] u (x, •) with the 
function u. By (6.1) we have that u t {x 1 t)dt = —(B(x,u)) x in the sense of measures, 



i.e. 



/ / 4>{x, t)u t (x, t)dt = - <t>{x, t) d (B(x, 



u)) x . 



Since the jumps of r] u (-,u)) are contained in J K , reasoning as in the proof of Corol- 
lary 4.4, we have that {rj u {x 1 u))*u t {x,t)dt = —(r] u (x,u))*(B(x,u)) x in the sense of 
measures, i.e. 

// 4>{x,t){q u {x,u))*u t {x,t)dt = - / / (p(x,t) (r) u (x,u))* d(B(x,u)) x . 

J JR.x]0,T[ J JRx]0,T[ 

Hence 

/ / <p(x, t) d [(f](x, u)) t + (q(x, u)) x ] 

J JRx]0,T[ 

= -// (f)(x,t)(r] u (x,u))* d(B(x,u)) x + <j>(x, t) d (q(x, u)) x , 

J Jlx]0,T[ J JRx]0,T[ 



A chain rule formula in BV 



23 



so that it is enough to prove that (rj u (x, u))*B(x, u) x > (q(x,u)) x in the sense of 
measures, i.e. for every nonnegative function g C c (Rx]0,T[) 



/ / <j)(x,t)(rj u (x,u))* dB(x,u) x 
Jo Ur 



dt > 



(j)(x,t) dq(x,u) a 



dt. 



We use the chain rule formula (see Corollary 4.4) and condition (E3) to obtain 

I ■= I I 4>{x,t){r) u {x,u))* dB(x,u) x dt 
Jo Ur 

/ 4>(x,t)r] u (x,u)B u (x,u(x)) ■ dDu(x) dt 
Jr 

I y2(j>(x,t)(r) u (x,u))*[B(x + ,u(x + )) - B(x_,u(x_))j 

JO ^ T . V 7 



-X£j K 



dt. 



We remark that the last term vanishes by the Rankine-Hugoniot condition. Using 
(E4) we obtain 



/ > 



T r 



dt 



{x,t)q u (x,u(x)) ■ dDu{x) 
+ / [ ^4>{x,t){q(x + ,u(x + )) - q(x-,u(x-))j 



dt 



t r 



cl>(x,t)dq(x,u) a 



dt. 



This concludes the proof. □ 

We consider the partially adapted Kruzkov entropies introduced by Audusse and 
Perthame for discontinuous flux (see formula (1.3) in [4]). 

In addition to the assumptions on the function B stated in Theorem 4.1, we also 
assume that 

for every x € K, the map B(x 7 •) is a one to one function from R to M . (6.4) 

Given a 6 1, by assumption (6.4) there exists a unique function c a : D a — > M, 
defined on a (possibly empty) set D a C K, such that B(x 7 c a (xj) = a for every 
x e D a . 

Proposition 6.2. For every a e K such that c a is defined in R, let us define 
the adapted Kruzkov entropy 

■q( a \x,u) := \u - c a (x)\ 

and the corresponding flux 

q {a) (x,u) := (B(x,u) - a)(sgn(u - c a {x)))* . 
We assume that for every x € Jk and for every (u_, u + ) e A x we have 

{sgn(u- - c a {x-)))* = (sgn(u + - c a (x+)))* . (6.5) 

Then we have 
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is an entropy-flux pair; in particular, the entropy inequality 

d t \u - c a (x)\ + d x [(B(x,u) - a)(sgn(u - c a (x)))*] < (6.6) 

holds in the sense of distributions; 
• (b) (6.2) holds for every entropy-flux pair (r],q) if and only if (6.6) holds for 
every a as above. 

Proof. For every i€l and u ^ c a {x) one has q^\x, u) = (sgn(u— c a (x))* B u (x, u). 
Then, since {r]C\x, it))* = (sgn(u— c a (x)))*, we obtain 7]^\x, u)B u {x, u) = q£*\x,u). 
Moreover, for every x € R and every (ii-,u + ) G A x satisfying (6.5) we have that 

q( a) (x + ,u + ) - q (a) (x-,u_) = (sgn(u+ - c a (x)))* (B(x+, u+) - a) 

— (sgn(u_ — c a (x)))* (B{x-,U-) — a) = . 

In order to prove (b), let u be a bounded BV solution to (6.1). If u satisfies (6.2) 
for every entropy-flux pair (r],q), then from (a) it satisfies also (6.6) for every a. 

Conversely, assume now that u satisfies also (6.6) for every a. Let (rj, q) be an 
entropy-flux pair, and let 4>(x, t) be a non-negative test function. We have to prove 
that (6.3) holds. 

Assume that \u(x, t) \ < M for every (x,t), supp0 c]a, 6[x]0,T[, and \B(x,u)\ < C 
for every (x, u) G]a, b[x] — M,M[. Let us fix a positive integer number N, and for 
every x £ [a,b] define 

I(x) :={ieZ: \i\ < N, a? := iC/N £ RangeB(x, •)}, 
m(x) := min/(a;), n(x) := max/(x), 
cf (x) :— c a N (x), i — m(x), . . . , n(x). 

We are going to approximate r? (and so q) by an entropy 77^ of the form 

n(x) — 1 

11 N {x,u):=a N {x)+b N {x)u+ £ b? (x)\u - cf (x)\ 

i—m(x)-\-l 
n(x) — 1 

= a N (x)+b N (x)u+ b?(x)r]^(x,u), 

i=m(x) + l 

where bf(x) > for every i and x. Indeed, if we define 



(6.7) 



5 i^ x ) - ^7 N ^(x) -c^(x) " " ' XER ' '' = " >{J ' ] " (- ' ] ' 1 



and 



T](x,cf +1 (x)) - r]{ x,cf{x)) 



b N (x) := 

, N( , 6^(x)-5^ 1 (x) . \ 

bi (x) := , i = m(x) + 1, . . . , n(x) - 1, 

n(x) — 1 

a N (x) := n(x, c N m(x) {x)) - b N ' (x)<% (x) (x) - £ b?{x)[c?{x) - <% (x) {x)], 

i=m(x)-\-l 
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then bf(x) > and rj N (x,-) is a convex piecewise afRne function coinciding with 
T](x, •) in the points Uj = (x), i — m(x), . . . , n(x). 
The flux associated to rj N is the function 

n{x) — 1 

q N {x,u) :=b N {x)B(x,u)+ ^ 6f {x)q«\x, u) . (6.8) 

i— m(x) + l 

For every JVeNwe have that 

/ [ 4>(x,t)[(7] N (x,u(x.t))) t + {q N (x,u(x,t))) x ]dx dt 

J Jrx[0,T] 

= cp(x,t)b N (x)[u t (x,t) + (B(x,u(x,t))) x ]dx dt 

J ilx [0,T] 

+ V // 0(x,t)^(a;)[(» 7 ( a « W )(x,«(x,t)))t + («^>(a;,^a;,t))) I ]didt, 

l=m (z)+l j ^xIO.T] 

(6.9) 

where we recall that 

?7 lai >(x, u) = \u — c a N{x)\ 

and 

q (a ?\x, u) := (B(a;, u) - c Q « (x)) (sgn(u - c a N (x))* . 

We recall that, given a non-negative measure fi (i.e. ^(</>) > for every test 
function (f> > 0), if we define a measure 



Mf>(0) : = / / (j)(x,t)b(x,t)dfi(x,t) , 

J JRx [0,T] 

where & is a non-negative Borel function, then is also a non-negative measure (see 
[15, Ch. 7]) . Hence from (6.1), (6.6) and (6.9) and the fact that the functions bf are 
non-negative, we have 



// 



(f>(x, t)[(r) N (x, u(x.t))) t + (q N (x, u(x, t))) x ] dxdt<0. 

T] 



The sequences of functions {ti n )n 1 {q N ) n are uniformly bounded on [a, b] x [— M, M], 
and converge pointwisc to rj and q respectively, hence we conclude that (6.3) holds. □ 
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